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Введение. Статья посвящена исследованиям термофлуктуаци- 
онных процессов в изоляционных материалах в соответствии 
с теорией теплопроводности для решения задач диагностики и 
прогнозирования остаточного ресурса изоляционных матери- 
алов на основе цифрового регистратора, а также неразруша- 
ющего температурного метода. 

Целью работы является совершенствование методов неразру- 
шающей диагностики, а именно — разработка автоматизиро- 
ванной системы контроля состояния изоляции и расчетно- 
экспериментальное исследование. 

Материалы и методы. Предложены математические модели, 
описывающие послойное распределение температуры кабель- 
ной линии в соответствии с теорией теплопроводности, ис- 
пользующие дифференциальное уравнение Фурье. Создан 
обобщенный алгоритм работы регистратора контроля пара- 
метров силовых кабельных линий (СКЛ), реализующий метод 
неразрушающего контроля термофлуктуационных процессов 
в изоляционных материалах СКЛ. Проведен сравнительный 
анализ экспериментальных и расчетных характеристик рас- 
пределений температуры. При этом исследовались различные 
нагрузочные режимы работы и функции изменения тока кабе- 
ЛЯ. 

Результаты исследования. Разработаны математические мо- 
дели и программное обеспечение для численного моделиро- 
вания температурного поля в сечении кабеля в соответствии с 
теорией теплопроводности, учитывающие физические свой- 
ства материалов и геометрические размеры элементов кабеля. 
Проведен сравнительный анализ экспериментальных и рас- 
четных характеристик распределений температуры. Научно 
подтверждена адекватность разработанной упрощенной ма- 
тематической модели определения температуры наиболее 
нагретой точки изоляции жилы кабеля на основе измеренных 
значений температуры поверхности силового кабеля и темпе- 
ратуры воздуха при различных изменениях действующего 
значения тока кабеля. Разработан и обоснован метод исследо- 
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вания термофлуктуационных процессов на основе использо- 
вания послойно расположенных в СКЛ температурных датчи- 
ков, позволяющий объединить в одном измерительном сред- 
стве два метода контроля: прогнозирование развивающегося 
дефекта изоляции и неразрушающий контроль термофлуктуа- 
ционных процессов силового кабеля. Предложенная матема- 
тическая модель может быть использована в качестве базовой 
при расчете тепловых процессов силовых кабелей в режиме 
реального времени, поскольку ее адекватность подтверждена 
экспериментальными исследованиями. 

Обсуждение и заключения. Полученные результаты могут 
быть использованы при разработке теории, методов диагно- 
стики и прогнозирования состояния изоляционных материа- 
лов в сложных распределенных системах при различных ре- 


жимах работы. 
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Введение. Увеличение перетока электроэнергии в энергосистемах приводит к тому, что СКЛ используют на 


пределе их физических возможностей, следствием чего 


является пробой электрической изоляции [1]. Таким образом, 


вопросы обеспечения энергетической безопасности, а также повышения эффективности передачи электроэнергии вы- 


ходят на первый план. 


Важной характеристикой СКЛ является нагрузочная способность, определяемая тепловым режимом. По ней 
рассчитывается площадь поперечного сечения и пропускная способность СКЛ [2]. Ресурс изоляционных материалов 
зависит от таких параметров, как температура, ток, влажность [3]. СКЛ с изоляцией из сшитого полиэтилена (СПЭ) 
обладает высокими термоизолирующими свойствами. Поэтому пропускная способность и надежность в значительной 
степени зависят от температурного режима работы. Для решения указанной проблемы используется мониторинг тем- 
пературы кабеля во время работы. Для исследования температурного поля в сечении СКЛ используется цифровой 


регистратор (рис. 1). 
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Рис. 1. Схема цифрового регистратора для проведения экспериментальных исследований 


Вт. 1. ДХлаэтат ор @еиИа/ гесогаег юг сопаисйиз ехретйтеп 5ше5 


Дубяго М. Н. и др. Метод исследования термофлуктуационных процессов в задачах диагностики 





Современные средства мониторинга распределенной по длине линии температуры не имеют подходящих для 
этого моделей и алгоритмов расчета тепловых процессов в изоляции кабеля [4]. 

Представленное исследование термофлуктуационных процессов основывается на методе регистрации темпе- 
ратуры внутри изоляции. Распределенные оптоволоконные системы с установленными в пределах поперечного сече- 
ния СКЛ волоконнооптическими датчиками измеряют внутренние температуры для определения точек повышенного 
нагрева. Для исследования законов передачи температуры в лабораторных условиях решено применить послойное 
размещение датчиков (рис. 1). 

При появлении частичного разряда (ЧР) в изоляционном материале рассеивается энергия, основная часть ко- 
торой затрачивается на деструкцию изоляционного материала. Поэтому в месте действия ЧР осуществляется нагрев 
изоляции, что и является объектом измерения [5, 6]. Учитывая общие требования к системам измерения термофлукту- 
ационных процессов и проблематику измерения, следует отметить, что существующие цифровые регистраторы не поз- 
воляют одновременно реализовать необходимый набор диагностических функций. Поэтому была использована система 
цифровой регистрации термофлуктуационных процессов, разработанная на кафедре электротехники и мехатроники 
(ЭимМ) Южного федерального университета (ЮФУ) [5]. 


Разработка алгоритма и системы мониторинга тепловых процессов СКЛ. Разработанная система измере- 
ния термофлуктуационных процессов сочетает в себе все достоинства современных цифровых измерительных систем, 
а именно: 

— мониторинг тепловых процессов с помощью системы 5САЛРА ‚ реализованной в среде ГаБИТЕЙ’; 
— мониторинг в режиме реального времени послойного измерения температуры; 
— регулирование и измерение влажности и т. д. (рис. 1 [5]). 

Исследование повышенного перегрева термофлуктуационным методом предусматривает использование в ла- 
бораторных экспериментах температурных датчиков, расположенных послойно в СКЛ (см. рис. 1), для создания мо- 
делей и алгоритмов расчета тепловых процессов в изоляции. 

Таким образом, разработанная схема автоматизированной системы контроля состояния изоляции 5МасСТС: 

— позволяет независимо или совместно использовать два метода контроля (прогнозирования развивающегося дефек- 
та изоляции и неразрушающего контроля термофлуктуационных процессов); 

— имеет набор диагностических функций по контролю режимов работы СКЛ (номинальной нагрузки, холостого хода, 
короткого замыкания) по результатам полученного профиля температур, снятых датчиками [5, 6]. 

Алгоритм работы программного обеспечения цифрового регистратора контроля параметров СКЛ представлен 
на рис. 2. 
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Рис. 2. Алгоритм функционирования регистратора контроля параметров СКЛ 
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Таким образом, обобщенный алгоритм: 
— во-первых, позволяет управлять процессом диагностики параметров СКЛ; 
—щ во-вторых, эффективно реализует метод контроля термофлуктуационных процессов в изоляционных материалах 
СКЛ посредством цифрового регистратора. 


Моделирование термопроцессов изоляции СКЛ с использованием пакета Гар /ТЕИ/. Для исследования тер- 
мофлуктуационных процессов в СКЛ создан виртуальный инструмент ГаБИЕЙ,, который: 
— позволяет взаимодействовать с измерительными датчиками и управляющей аппаратурой; 
— осуществляет сбор, обработку, отображение информации и результатов расчета; 
— моделирует как отдельные объекты, так и цифровую систему 5МаСТС в целом. 
Вид панели виртуального инструмента представлен на рис. 3. 
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Рис. 3. Панель виртуального инструмента ГаБ/ТЕ 
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Стенд для исследования термофлуктуационных процессов представляет собой компьютерную модель 
ГабБ ТЕТ? на рабочем столе персонального компьютера и содержит модели виртуального инструмента: 

— данные и графики температур; 
— данные и зависимости влажности; 
— анализ и средние данные температур; 
— временную зависимость тока [6]. 

Программное обеспечение позволяет рассчитывать значения тока, напряжения, температуры жилы кабеля, 
влажности, тока нулевой последовательности (тнп) с интервалом времени 1 с. 

Таким образом, создана автоматизированная система контроля характеристик изоляционных материалов 
5МастТС, позволяющая контролировать термофлуктуационные процессы СКЛ напряжением 10 кВ. При этом исполь- 
зуется метод неразрушающего контроля термофлуктуационных процессов в изоляционных материалах СКЛ. Специаль- 
ное программное обеспечение снимает информацию о работе $МаСТС. 


Мониторинг тепловых условий кабеля и окружающей среды. Для исследования тепловых режимов СКЛ 
использовались математические и имитационные модели силового кабеля с СПЭ [6, 7, 8, 9]. В этом случае кабель и 
окружающая среда разбиваются на зоны, представленные в виде электрической схемы замещения (рис. 4). 
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Рис. 4. Сечение кабеля (а) и его эквивалентная тепловая схема замещения (5) 
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Здесь 0,,, 9,, Он, 0. о, 9.с. — соответственно температуры на жиле, экране, основной изоляции, защитной обо- 
лочке, окружающей среде; С», С., Сои, Со, Сос — теплоемкости соответственно жилы, экрана, основной изоляции, 
защитной оболочки, окружающей среды; К», К., Во и, В.о— тепловые сопротивления соответственно жилы, экрана, 
основной изоляции, защитной оболочки; О„, О,„и— источники тока, замещающие температурный напор, создаваемый 
током, протекающим по жиле и экрану. 

На рис. 4 показаны слои и точки, в которых проводится контроль тепловых параметров. Исследование рас- 
пределения температуры кабельной линии выполнялось в соответствии с теорией теплопроводности, при этом ис- 
пользовалось дифференциальное уравнение Фурье [10]. 

Физические свойства материалов и геометрические размеры элементов кабеля АПвПу г-1х30/25-10 [5, 6] све- 
дены в табл. 1 [11]. 

Таблица 1 


ТаЫе 1 
Параметры кабеля 


СаЫе рагатеегз 


М 

Параметр 
Удельная теплопроводность, Вт/(м`К) 209,3 0,024 
о | 2% | 8100 





00 
Удельная теплоемкость, Дж/(кг-К) 1900 1005 


Профиль температур в сечении силового кабеля (рис. 4) определялся системой уравнений в соответствии с теорией 





теплопроводности. 

0 + (= - *—)/2пих, (= Вос) 
(и) + вт (=) (+ =)/2тжье, вы», В < Ве) 
О(, о.) + (°.) (+1) /2ийьи. „№ < В < йо) 


9(^) = 


1% 8512 п 2 2 


э п? 
9(#,) + м () 2ПАь 4^у. , (Йон. < А 1.) 


2 
ВСР) + ва ("9 = /2пА и, (к 5 В 5 йаи) 
9(н») + 52 (#2 — #2) /42жуж, (0 < В < Аь). 


Здесь Лос, Лои,Ль, Аж — Удельные теплопроводности окружающей среды, полиэтилена, экрана и жилы соот- 


ветственно; [„,[, — токи на жиле и экране; х — коэффициент теплопроводности; 6х = [«/5ж; 6. = 5/55; 5ю95— 


площади поперечного сечения жилы и экрана; у, У. — электропроводности жилы и экран. 
В соответствии с системой уравнений (1) в МаН.аб строится профиль температур в сечении кабеля и окружа- 
ющей среды [6, 12]. 


Сравнительный анализ экспериментальных и расчетных характеристик распределений температуры. 
Для доказательства адекватности предлагаемой упрощенной математической модели СКЛ авторы провели исследова- 
ния, в ходе которых через кабель пропускался переменный ток, величина которого изменялась в зависимости от ре- 
жима работы (рис. 5). 
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Рис. 5. Экспериментальная диаграмма тока в кабеле 
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В ряде работ [1, 6, 13] представлены модели, с помощью которых рассчитывается распределение температуры 
по слоям, а также зависимость, позволяющая определять температуру жилы (6х. °С). 

Предлагается упрощенная аналитическая зависимость (2) определения температуры наиболее нагретой точки 
изоляции (жилы кабеля), полученная из системы (1). Температура жилы кабеля вычисляется с учетом температуры 
поверхности кабеля, окружающей среды и тока жилы: 


п-12.р2о’т-Тк-Кр|1+а-(@доп-@о.с.(6)| 


0» (Е) = 9. (о + РЕ бен (2) 

В опытах использовался силовой кабель АПв Пу г-1х30/25-10. 

Здесь 9, — расчетная температура жилы кабеля, °С; 0. „.=23 °С — измеренная температура поверхности ка- 
беля и оболочка); и = 1 — число жил кабеля; Г — максимальный ток кабеля при проведении измерений, А; 
р»о= 2,8-10*Ом-м — удельное электрическое сопротивление жилы кабеля при 20 °С; Т, = 0,0028 °С, м/Вт — сумма 
термических сопротивлений изоляции и защитных покровов кабеля; и — экспериментальный множитель; К = 1,02 — 
поправочный коэффициент для приведения электрического сопротивления к расчетной температуре; а = 0,0043081 °С 
— температурный коэффициент сопротивления материала жилы [7]; 9 = 90 °С — длительно допустимая температу- 
ра изоляции жилы; 09... = 23 °С — измеренная температура окружающей среды; 5 = 0,0038465 — сечение жилы кабе- 
ля, М’. 

Итак, 
1.1?.2,8.10-8. 0,0028 . 1,02[1 + 0,004308 - (90—60... (Е))| 
0,„.(Е) = 9....(Е) + 38465.10“ 

Основное назначение уравнения (2) — рассчитать температуру жилы кабеля максимально просто. Уравнение 
учитывает лишь стационарные тепловые процессы в кабеле и справедливо при установившихся значениях измеряе- 
мых величин. Это характерно и для всех нормативных документов по расчету тепловых процессов в силовых кабелях: 
МЭК (ТЕС) 60287, ГОСТР МЭК 60287-1-1, ГОСТР МЭК 60287-1-2, ГОСТР МЭК 60287-1-3, ГОСТ Р МЭК 60287-2-1, 
ГОСТР МЭК 60287-2-2. 

При эксплуатации вполне естественно, что измеряемые величины изменяются во времени, таким образом, для 
их учета необходимо вести расчет нестационарных тепловых процессов. 

Эксперимент 1. Экспериментальная диаграмма ступенчатого изменения действующего значения тока кабеля 
приведена на рис. 5. Экспериментальные и расчетные временные зависимости процесса нагрева кабеля, выражение 
(2), построенные при ступенчатом изменении тока, приведены на рис. 6. 


[ИАМИАИААИ 


Температура, °С 


Рис. 6. Экспериментальные и расчетные временные зависимости процесса нагрева кабеля 
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В результате экспериментальных исследований получены показатели качества предложенной математической 
модели тепловых процессов силового кабеля: 
— относительная погрешность модели не превышает 11 %; 
— среднеквадратическое отклонение экспериментальных данных от расчетных по модели не превышает 4,5 °С. 

Не учтено влияние теплообмена через торцы исследуемого кабеля в исходной модели. Во время проведения 
эксперимента торцы жилы не были закрыты теплоизолирующей тканью, и происходил теплообмен с держателями. 
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Из-за малой длины образца этот теплообмен оказал влияние на результат. Приведенные диаграммы демонстрируют 
корреляцию экспериментальных и расчетных кривых температуры, что говорит о верном выборе математической мо- 
дели [2,4]. 

Эксперимент 2. Для проверки адекватности разработанной математической модели на основе мощного пони- 
жающего трансформатора с выходным током до 640 А была создана измерительная система (рис. 1), позволяющая 
исследовать тепловые процессы в одножильных кабелях. Датчики температуры послойно размещались внутри и на 
поверхности СКЛ, а также на некотором удалении от него. Значение тока определялось посредством измерительного 
трансформатора. Все измерения проводились с помощью цифрового регистратора, подключенного к компьютеру. В 
качестве образца использовались отрезки заведомо исправных кабелей с изоляцией из сшитого полиэтилена АПв 


Пу г-1х30/25-10 длиной 0,85 м. Экспериментальная диаграмма действующего значения тока приведена на рис. 7. 
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Рис. 7. Экспериментальная зависимость для действующего тока СКЛ 
Ето. 7. Ехретттета 4ерепаепсе юг РСЕ утх5 сиггепй 


На рис. 8 приведены экспериментальные и расчетные временные диаграммы температуры кабеля АПв Пу г- 
1х30/25-10 по опыту, представленному в рис. 7. 


Эспериментальня 


Расчетная 


Температура, °С 


ьс 
1494 2979 4482 7179 8662 10156 





Рис.8. Временные диаграммы температуры на жиле кабеля АПв Пу г-1х30/25-10 
Е 1.6. Тетрегаииге Иите @аэтат$ а са соте АРу Ри 8-1х30 / 25-10 


Расчет показателей качества тепловой модели СКЛ. Введем показатели качества тепловой модели для 
анализа и оценки погрешностей расчета температуры по рассматриваемым тепловым моделям [14]. 

Пусть (Д0); — абсолютная погрешность расчета температуры жилы на интервале расчета (или измерения) /. 
То есть 
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(46); = (6), — (6, 
где / — номер интервала расчета. 
В качестве показателя качества будем рассматривать величину 
ДА@млкс = тах!|(Д6) |}. 
Погрешность ДАд9мдкс характеризует диапазон точности тепловой модели, то есть границы, в которых могут 


находиться значения величин (Д6,,),. 


Рассмотрим также средние значения векторов (Дб); и (460,). как показатель качества тепловой модели: 
7 


п 
1=1 |(д6,) В 
р п ; 
— Х/=1[(А6.); | 
п 
где п — количество интервалов на расчетном промежутке времени. 
Величины ДО, И 6, характеризуют ошибки расчета температуры жилы в статике, то есть по ПОСТОЯННОЙ со- 


ставляющей. Для анализа качества тепловой модели в динамике будем рассчитывать коэффициент корреляции: 


1650) — 6, - У [6( — 8.2 





где 
в — 1 0р0) Аб — 2-1 050) 
р п о п | 
Обсуждение и заключения. Результаты расчета показателей качества математической модели по экспери- 
менту 1, полученные путем сравнения измеренной и расчетной температуры кабеля, приведены в табл. 2. 
Таблица 2 
ТаЫе 2 


Показатели качества тепловой модели при ступенчатом изменении тока 


Трегта! тоае! диаПё? рагатеегу ипаег “ерил5е сиггепт! уатаНоп 


Значения погрешности по формуле 
Показатели точности расчета - - 
стандартной предлагаемой 
Диапазон точности Адмдкс 22,419 7,657 


Показатель качества (среднее значение векторов) 
Коэффициент корреляции 
Максимальная относительная погрешность, % ООО: _ЗИОООООЯ СОООООООИАИЮ С СЗЗОООООИАИА 
Среднеквадратическое отклонение 


Полученные высокие значения коэффициента корреляции свидетельствуют о соответствии математической 





модели реальным тепловым процессам. Максимальная абсолютная погрешность при сравнении экспериментальных 
данных и расчетов по исходной формуле (рис. 6) составила 22,4 °С. Максимальная абсолютная погрешность при срав- 
нении экспериментальных данных и расчетов по предлагаемой формуле (рис. 6) составила 7,6 °С. 
Результаты расчета показателей качества математической модели по эксперименту 2, полученные путем 
сравнения измеренной и расчетной температуры кабеля, приведены в табл. 3. 
Таблица 3 
ТаЫе 3 


Значения показателей качества тепловой модели 


Трегта! тоае! диа р рагатеегу уа[ие5 


БИр://уезшК.Аоп$а.га 


Коэффициент корреляции 0,999 
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Максимальная абсолютная погрешность расчета температуры по отношению к непосредственно измеренной 
температуре не превысила 2,9 °С. Высокие значения коэффициента корреляции свидетельствуют о соответствии ма- 
тематической модели реальным тепловым процессам. 

Предложенная аналитическая зависимость может применяться в качестве базовой при вычислении тепловых 
процессов СКЛ в режиме реального времени, так как ее адекватность подкреплена опытными данными. 


Выводы. Исследования показали, что разработанные расчетные соотношения даже в условиях принятых 
упрощений достаточно корректно описывают тепловые процессы в кабеле. Применение рассмотренных устройств 
повысит эксплуатационную надежность систем электроснабжения. Основная область применения разработанной ма- 
тематической модели — диагностика и прогнозирование ресурса изоляции силового кабеля. Цифровое аппаратное 
устройство способно измерять температуры поверхности кабеля и окружающей среды, а затем в режиме реального 
времени рассчитывать внутренние температуры и решать задачи раннего выявления повреждений. 

Работа выполнена по гранту № ВнГр-07/2017-15 «Разработка теоретических основ и методов построения ин- 
теллектуальных многосвязных систем управления процессами производства, транспортировки, распределения и по- 
требления энергии». 
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